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Abstract 

A generalization of the parton-like formula is used for the first time to find the 
differential distributions in the inclusive semileptonic weak decays of the B meson. 
The main features of this new approach are the treatment of the 6-quark as an 
on-mass-shell particle and the inclusion of the effects arising from the 6-quark 
transverse motion in the S-meson. Using the 6-quark light-front (LF) distribution 
function related to the equal time momentum wave function taken from the ACM 
model we compute the electron energy spectra and the total semileptonic widths of 
the B meson. We find an impressive agreement between the electron energy spectra 
calculated in the LF approach and the ones obtained in the ACM model, provided 
the b quark mass is identified with the average of the floating b quark mass in the 
ACM model. In spite of the simplicity of the model we obtain a fair good description 
of the CLEO data for \Vcb\ = 0.042. 
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1. Introduction. Studying of the leptonic spectra in the weak semileptonic B —>■ X(.{Xu)C-i'e, 
decays^ is important to extract \Vcb\ and \Vub\- The latter plays an important role in the 
determination of the unitarity triangle. Since the 6-quark is heavy compared to the QCD 
scale, the inclusive semileptonic B decays can be treated with the help of an operator 
product expansion (OPE) combined with the heavy quark expansion (HQE) [|I|. The 
result (away from the endpoint of the spectrum) is that the inclusive differential decay 
width dr/dE may be expanded in A/nib, where A is a QCD related scale of order 500 
MeV, and nib is the mass of the heavy quark. The leading term (zeroth order in A/mb) 
is the free quark decay spectrum, the subleading term vanishes, and the subsubleading 
term involves parameters from the heavy quark theory, but should be rather small, as it 
is of order (A/m^)^. However, near the end point the l/rub expansion has to be replaced 
by an expansion in twist. To describe this region one has to introduce a so-called 'shape 
function', which in principle introduces a large hadronic uncertainty. This is quite analo- 
gous to what happens for the structure function in deep-inelastic scattering in the region 
where the Bjorken variable — 1. A model independent determination of the shape 
function is not available at the present time, therefore a certain model dependence in this 
region seems to be unavoidable, unless lattice data become reasonably precise. 

As to phenomenological analyzes of the photon spectra up to now they have been 
solely based on the ACM model 0, 0]. The various light-front (LF) approaches to 
consideration of the inclusive semileptonic transitions were suggested in Refs. In 
Refs. [0, the Infinite Momentum Frame prescription ph = C,Pb, and, correspondingly, 
the floating b quark mass ml{C,) = ^^M^ have been used. The transverse b quark momenta 
were consequently neglected. In Ref. the 6-quark was considered as an on-mass-shell 
particle with the definite mass rub and the effects arising from the 6-quark transverse 
motion in the B-meson were included. The corresponding ansatz of Ref. reduces to 
a specific choice of the primordial LF distribution function \4'{^yP'±)\'^ y which represents 
the probability to find the b quark carrying a LF fraction ^ and a transverse momentum 
squared = |p±p. As a result, a new parton-like formula for the inclusive semileptonic 
b ^ c,u width has been derived []^, which is similar to the one obtained by Bjorken et 
al. in case of infinitely heavy b and c quarks. 

In this paper, we use the techniques developed in Ref. [[7| to evaluate the non- 
perturbative corrections to the lepton spectrum in the inclusive B — > Xctvi decays. We 
strive to implement the binding and the S-meson wave function effects on the lepton 
energy spectrum. The main purpose of the paper is to confront the lepton spectra and 
the total semileptonic decay widths calculated in the LF and ACM approaches. We find 
that the discrepancy between the two is very small numerically. 

2. ACM versus LF. The decay spectrum in the ACM model is determined by the kine- 
matics constrains on the b quark. It incorporates some of the corrections related to the 
fact that the decaying b quark is not free, but in a bound state. The model treats the B 
meson with the mass ms as consisting of the heavy b quark plus a spectator with fixed 
mass m^p; the latter usually represents a fit parameter. The spectator quark has a mo- 
mentum distribution 0(|p|) where p is its three-dimensional momentum. The momentum 

^Throughout this paper we use the charge conjugated notations. 
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distribution is usually taken to be Gaussian: 




which is normalized so that the integral over all momenta of 0(|p|)p^ is 1. The energy- 
moment uii]_cons_ervat ion in the B meson vertex implies that the b quark energy Ef, = 
niB — i/p^ + ^Ip, where is the mass of the B meson; thus the h quark cannot possess 



a definite mass. Instead, one obtains a "floating" h quark mass {mlY = ra\ + m^^ — 
2mB^p2 + mlp , which depends on |p|. The lepton spectrum is first obtained from the 

spectrum dTf* (m/, E') / dE' of the h quark of invariant mass m{ (in the h quark rest frame) 
and then boosting back to the rest frame of the B meson and averaging over the weight 
function 0(|p|) (further details can be found in the original work P). 

The LF approach differs from that of the ACM model in two respects. First, similar to 
the ACM model the LF quark model treats the beauty meson B as consisting of the heavy 
h quark plus a spectator quark. Both quarks have fixed masses, and m^p, though. This 
is at variance with the ACM model, that has been introduced in order to avoid the notion 
of the heavy quark mass at all. Secondly, the calculation of the distribution over lepton 
energy in the LF approach does not requires any boosting procedure but instead is based 
on the standard Lorentz-invariant kinematical analysis (see e.g. Ref. [^]), which we now 
briefly discuss. 

There are three independent kinematical variables in the inclusive phenomenology, 
for which we choose the lepton energy -E^, g^, where q = Pe + Pug, and the invariant mass 
~ {Pb—qY of a hadronic state. Introducing the dimensionless variables y = 2E^/mB, 
t = /m\, and s = Mx/m^, the differential decay rate for semileptonic B decay can be 
written as 

'^^^ - ^^^Vc.P f dt f ds 



dv QAtc^ 

^ so 

X \tW^ + i[y(l + t - s) - _ t\W2 + t[ ^^\~' - 2/]W^3 + ...}, (2) 



where the structure functions Wi = Wi{s,t) appear in the decomposition of the hadronic 
tensor Wa/B in Lorentz covariants 0. The ellipsis in (Q) denote the terms proportional 
to the lepton mass squared. The kinematical limits of integration can be found from 
equation 

- < 1- (3) 



1 - y y 

, 2 



They are given by < ?/ < i/max = 1 — ("^x^*'^V''^b) 5 where rrix^^^^ is the minimal mass 
of the hadronic state X^, s^ax = ^ + t-iy + t/y), and tmax = y ((1 - (1 - ymax)/i^ - y))- 



3. Light-cone distribution functions. In a parton model we treat inclusive semileptonic 
B Xcii'i decay in a direct analogy to deep-inelastic scattering. Specifically, we assume 
that the sum over all possible charm final states can be modeled by the decay width 
of an on-shell b quark into on-shell c-quark weighted with the 6-quark distribution. 
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Following the above assumption, the hadronic tensor Wcb is written as 



where 



Ur 



2 -^^ ui-u --u-'/3"'C) 

spins 



(4) 



(5) 



with Oa = 7a(l — 75)- The factor 1/^ in Eq. (^) comes from the normalization of the B 
meson vertex [0. 

Eq. (^) amounts to averaging the perturbative decay distribution over motion of heavy 
quark governed by the distribution function f{x,p\) = Iii this respect our 

approach is similar to the parton model in deep inelastic scattering, although it is not 
really a parton model in its standard definition. The normalization condition reads 



TT 



1 



(6) 



The function 6{ec) where Ec is the c-quark energy is inserted in Eq. for consistency 
with the use of valence LP wave function to calculate the 6-quark distribution in the 
i?-meson. 

Since we do not have an explicit representation for the B-meson Fock expansion in 
QCD, we shall proceed by making an ansatz for 'ip{C,,p'±)- This is model dependent 
enterprise but has its close equivalent in studies of electron spectra using the ACM model. 
We choose the momentum space structure of an equal time (ET) wave function 0(|p|) 
as in Eq. (|I|). We convert from ET to LF momenta by leaving the transverse momenta 
unchanged and letting 



1 



Pi 



(Pf -Pi ) = o(Pi 



pf ' 



(7) 



for both the 6-quark [i = b) and the quark-spectator [i = sp). The longitudinal LF 
momentum fractions are defined as ^sp = Ptp/^B^ = pt I with ^b + isp = 1- In 
the i?-meson rest frame Pg = tub- Then for the distribution function \il'{^,p'j_)\'^ 
normalized according to (j^) one obtains 



6) 



l^(e,pi)P 



^^exp 



pI + pV 

P) , 



dpz 



where 




p]_ + mf 



sp 



j and 


dpz 







TUB + 



pi + ^% 



iYrriB^ 



(9) 



^On the same footing one can consider |p| in (|l|) as the relative momentum between heavy and Ught 
quarks. In this case it is more convenient to use the quark-antiquark rest frame instead of the B- 
meson rest frame. Recall that in the LF formalism these two frames are different. Then the longitudinal 
LF momentum fractions are defined as S^sp — ptp/'^'^0 7 Cfc — pt/-^^o, where the free mass Mq is 

p2 with 



Mo = ^Jml + p2 + sjm% + p2 with p^ = p\+pl, = 

explicit form of \dpz/d^\ is given e.g. by Eq. (10) of Ref. 
approaches yield identical results for the electron spectra 



i)Mo — {ml — mlp)/2Mo. In this case the 
. We have checked that numerically both 
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The calculation of the structure functions Wj(t, s) in the LF parton approximation (^) 
is straightforward. The result is 

W,{t,s) = I w,is,t,^,pl)S[ip, - qf - ml] ^'^ (eJ^e^V, (10) 

where Wi{s, t, ^, p\) are the structure functions for the free quark decay. For further details 
see Appendix of 0. Eq. (|1^) differs from the corresponding expressions of Refs. [^] and 
by the non-trivial dependence on p\ which enters both and argument of 

the 5-function. For further details see 

4- Results. Having specified the non-perturbative aspects of our calculations, we proceed 
to present numerical results for the lepton spectrum in the decay B — > Xceu^. Our main 
computation refers to the case m^p = 0.15 GeV, rric = 1.5 GeV as chosen in Ref. [O . 



The choice of mf, in our approach deserves some comments. In the ACM model, it 
was shown [|ri[|, |[T^ (see also |]T3[) that the corrections to first order in l/rrib both to 
the inclusive semileptonic width and to the regular part of the lepton spectrum can be 
absorbed into the definition of the quark mass: m^*"*^ =< ml >, where < m( > is the 
value of the floating mass, averaged over the distribution 0(|p|). 

The choice of mi, in the LF approach was first addressed in the context of the LF model 
for b ^ s'y transitions in Ref. |T^. Using the scaling feature of the photon spectrum in the 



LF model, it was suggested that m^^ can be defined from the requirement of the vanishing 
of the first moment of the distribution function. This condition coincides with that used 
in HQE to define the pole mass of the 6-quark. In this way one avoids an otherwise large 
(and model dependent) correction of order l/mi, but at expense of introducing the shift 
in the constituent quark mass which largely compensates the bound state effects. It has 
been also demonstrated that the values of m^^ found by this procedure agree well with 
the average values < m( > in the ACM model. The photon energy spectra calculated in 
the LF approach were found to agree well with the ones obtained in the ACM model. 

Accepting the identification m^^ = m^'~^^^ we want now to check whether this result 
holds for the description of the other channels, like 6 — c. We find again a good agreement 
between the LF and ACM results but now for the semileptonic 6 — c decays. The results 
of our computations of electron spectra and the semileptonic widths are reported in Fig.l 
and Table 1 using |V"cfe| = 0.04. The different curves in Fig.l correspond to the different 
values oipf. For each case we show separately the inclusive differential semileptonic decay 
widths for the LF and ACM models and the free quark decay. The quark decay spectra 
vanish for y > {mb/mB){l — ml/ml), whereas the physical endpoint is Umax = 1 — "^d/^^b, 
where mo is the mass of the D meson. In the LF approach the endpoint for the electron 
spectrum is in fact not Umax but y^ax = 1 ~ ml/m%. This is the direct consequence of 
the p\ integration in Eq. (|10D 0. Note that coincides with with accuracy 

~ msp/mB- For m^ ~ 1.5 GeV the difference between y^^ ~ V^x' ymax is of the 
order 10~^. Another possibility advocated in is to sum the electron spectra from the 
exclusive B ^ D,D* channels and from the inclusive B ^ Xc channels, where Xc is the 
hadronic state with the mass mx^ > mu*. Such the 'hybrid' approach will be considered 
elsewhere. 

In Table 1 for various values of pf, we give the corresponding values of the total 

4 



semileptonic width for the free quark with the mass m^ =< ml > and the B meson 



^Note that the expressions for Wi in Ref. Q miss an extra factor of 2. 
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semileptonic widths, calculated using the LF and ACM approaches, respectively. In the 
last two columns, we give the fractional deviation 6 = ATsl/^sl (i^ P^r cent) between 
the semileptonic widths determined in the LF and ACM models and that of the free quark. 
The agreement between the LF and ACM approaches for the electron spectra is excellent 
for small pf as is exhibited in Fig.l. A similar agreement also holds for integrated rates 
shown in Table 1. This agreement is seen to be breaking down at pj > 0.4 GeV, but even 
for Pf ~ 0.5 GeV the difference between the ACM and LF inclusive widths is still small 
and is of the order of a per cent level. 

Finally, we calculate the b —>■ c spectrum and compare it with the experimental data 
from the CLEO collaboration . This is a direct calculation of the spectrum and not a 



X fit. We briefly investigated the sensitivity of the electron spectra to other parameters 
of the models and found that the choice pf = 0.4 GeV, rric = 1.5 GeV rrisp = 0.15 GeV is 
quite acceptable. 

We display the results in Fig. 2, where the three theoretical curves are presented for the 
LF, ACM and free quark models. In these calculations we have implicitly included the 
0{as) perturbative corrections arising from gluon Bremsstrahllung and one-loop effects 



which modify an electron energy spectra at the partonic level (see e.g. [|T6[ and references 
therein). It is customary to define a correction function G{x) to the electron spectrum 
dr^^^ calculated in the tree approximation for the free quark decay through 

f^'^h- |^G(.)) , (11) 

ax ax \ Sir J 

where x = 2E/mb. The function G{x) contains the logarithmic singularities ~ ln^(l — x) 
which for rric = appear at the quark- level endpoints Xmax = 1- This singular behaviour 
at the end point is clearly a signal of the inadequacy of the perturbative expansion in this 
region. The problem is solved by taking into account the bound state effects 0. Since 
the radiative corrections must be convoluted with the distribution function the endpoints 
of the perturbative spectra are extended from the quark level to the hadron level and the 
logarithmic singularities are eliminated. In actual calculations we neglect the terms ~ p 
in G{x) and take this function from |16|, Eq. (4.10). 



The agreement with the experimental data is good. Using |V"c6| = 0.042 for the overall 
normalization^ we obtain for the semileptonic branching ratios 

BRlf = 10.16%, BRacm = 10.23% BRfree = 10.37%, (12) 
in agreement with the experimental finding |]T5| BRsl = 10.49 ± 0.17 ± 0.43%. 



4. Conclusions. We have applied a new LF formula ^ to calculate the partial electron 
spectra in the semileptonic B decays. Using the ET and LF 6-quark distribution functions 
related by a simple kinematical transformation we compared the LF and ACM models by 
computing the 6 — > c decays for m^-^ =< m{ >. A summary of our results is presented 
in Table 1 and Fig.l shows a good agreement between the results of the two models. We 
have also calculated (Fig. 2) the b —>■ c spectrum including the perturbative corrections 



and found an agreement with the experimental data from the CLEO collaboration |15 



A more detailed fit to the measured spectrum can impose constrains on the distribution 

"^This value agrees with the combined average \Vcb\ — (40.2 ± 1.9) x lO^'^ of Ref . p^. 
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function and the mass of the charm quark. Such the fit should also account for detector 
resolution. 

The same formulae can be also applied for nonleptonic B decay widths (corresponding 
to the underlying quark decays b cqiq2) thus making it possible to calculate the B 
lifetime. A preview of this calculation can be found in [jl8[. It would be interesting to 



check whether the effective values of the 6-quark mass can appear to be approximately 
the same for different quark channels and for different beauty hadrons. This work is in 
progress, and the results will be reported elsewhere. 



We thank Marco Battaglia and Pepe Salt for the discussion and Karen Ter-Martirosyan 
for his interest in this work. This work was supported in part RFBR grants Refs. 00-02- 
16363 and 00-15-96786. 



Table. For the values of p/ in column (1) we display the average value of the floating 
b-quark mass < m( > (both in units of GeV) in the second column and the total semilep- 
tonic width of the free b-quark (in units of ps~^) in the third column. In the forth and 
fifth columns we compare the total semileptonic widths, calculated in the ACM and LP 
approaches, respectively. In all cases m^p = 0.15 GeV and nic = 1.5 GeV and the radia- 
tive corrections are neglected. In the sixth and seventh columns, we give the fractional 
deviation in percent between the semileptonic widths determined in the LP and ACM 
models and that of the free quark. A momentum distribution of the b-quark is taken in 
the standard Gaussian form ([^) with the Permi momentum pf. \Vcb\ = 0.04. 
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Figure 1. The inclusive differential semileptonic decay widths dTsi/dEe in unites 
ps~^GeV~^ for the LF model (thick solid lines), ACM model (thin solid lines), and the 
free quark decays (dashed lines). The parameters are pf = 0.2 — 0.5 GeV, rngp = 0.15 
GeV, and = 1.5 GeV. I^,,! = 0.04. 
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Figure 2. The predicted electron energy spectrum compared witli tlie CLEO data |T5 . 
Tlie calculation uses pf = 0.4 GeV, mh = 4.8 GeV, rric = 1.5 GeV, and = 0.25 for 
the perturbative corrections. Thick solid line is the LF result, thin solid line is the ACM 
result, dashed line refers to the free quark decay. The spectra normalized to 10.16%, 
10.23%, and 10.36%, respectively. \Vcb\ = 0.042. 
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